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Abstract. The multiconfiguration Dirac-Hartree-Fock (MCDHF) model has been 
employed to calculate the transition rates between the nine lowest levels of radium. 
The dominant rates were then used to evaluate the radiative lifetimes. The decay 
of the metastable 7s7p 3 Po state through 2-photon E1M1 and hyperfine induced 
channels is also studied. 

1 Introduction 

Recent advances in trapping and spectroscopy of free, neutral atoms make it possible to extend 
the search for time reversal violation effects into the domain of radioactive elements [T] . In the 
last decade several heavy atoms were considered as candidates for experimental searches [2]. 
There arc at least two ongoing atomic trap experiments (in Kernfysisch Versneller Instituut 
in Netherlands [314] and in Argonne National Lab in the U.S. [5l6p whose aim is to detect 
the electric dipole moment of radium. The advantage of radium lies in octupole deforma- 
tions of nuclei in several isotopes [7J, simple electronic structure (ground state configuration 
[Kr]4d 10 4/ 14 5s 2 5p 6 5d 10 6s 2 6p 6 7s 2 yields the closed-shell singlet state ^So) as well as in coin- 
cidental proximity of two atomic states of opposite parity, 7s7p 3 P\ and 7s6d 3 D2, which are 
separated by a very small energy interval 5 cm -1 . 

The data on atomic spectrum of radium compiled in the tables of Moore [8] came from 
the experimental investigations of Rasmussen [9], with subsequent revisions by Russell [ID] : 
both go back to the 1930s. These data cover only 69 classified lines. The isotope shifts and 
hyperfine structures of radium were measured by the group of Wendt [11|12|13|14|15|16] in 
the 1980s. They have studied both atomic Ra I and ionic Ra II spectra, and obtained the 
isotope shifts, magnetic dipole hyperfine constants A and electric quadrupole constants B of 
the 7s7p 1 P\, 7s7p 3 P\, 7s7p 3 P 2 , and 7s7d 3 D 3 levels of neutral radium, as well as the magnetic 
dipole hyperfine constants A of the 7s 2 S\/2, 7p 2 Pi/ 2 , and 7p 2 P 3 / 2 levels of Ra + ion, together 
with the electric quadrupole constant B of the 7p 2 P 3 / 2 state of Ra + . Hyperfine structures of 
singly- ionised radium have been also the subject of several theoretical papers [17|18 19 20 2 1|22j . 
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The excitation energies of several states of neutral and singly-ionised radium (and barium) 
were later calculated in the framework the relativistic coupled-cluster theory by the group 
of Kaldor [23124] . More recently there have been three papers from our group in which we 
calculated the lifetimes of the 7s7p x Pi and 7s6d 3 D^ states [25] ; the hyperfine structure constants 
of all levels belonging to the two lowest excited-state configurations and the electric field 
gradients generated by the electronic cloud in the 7s7p 1 Pi, 7s7p 3 Pi, and 7s7p 3 P2 states, which 
in turn (combined with the measured values of the electric quadrupole constants B) yielded 
the nuclear electric quadrupole moment of radium-223 isotope |27j . The excitation energies and 
lifetimes of several states of radium (and barium) were calculated recently in the framework 
of the combined configuration-interaction and many-body perturbation theory by Dzuba and 
Ginges [2S]. 

In the present paper we calculated the transition probabilities between the states arising 
from the three lowest configurations of radium: 7s 2 ^So, 7s7p 3 Po,i,2i 7s7p l P\, 7s6d 3 Di,2,3, and 
7sQd 1 I?2, as well as the lifetimes of these states. The purpose of the present paper is fourfold. 
Firstly, we intended to extend the transition rate calculations on all levels arising from the three 
lowest configurations (i.e. 7s 2 , 7s7p, and 7s6e?). We included (1) transitions which contribute 
appreciably to lifetimes; (2) transitions involving 7s6d and 7s 7p 3 Pi levels because they are 
of interest in EDM experiments (the mixing induced by T-odd interactions is strongest between 
these two levels); (3) transitions which may be important for trapping; (4) transitions which are 
stronger than 1/s (a somewhat arbitrary threshold). Secondly, new comparison of our results 
became available when the two most recent papers [2815] appeared in print. These in turn 
permitted further tests of the newly developed [29] parallel version of the GRASP package [30] 
and calibration the theoretical model for the calculations of the spectroscopic properties of 
radium. Finally, we present a summary of the available theoretical data on the spectroscopic 
properties of eight lowest excited levels of radium, with the hope that they may be of help for 
the experimental groups, that are currently in the process of setting up the atomic traps for 
the search of permanent electric dipole moments. The Argonne group has already measured [5] 
the frequency and the rate of the 3 P\ — x So transition in the 22 jRa isotope, and determined the 
lifetime of the 3 P\ level. The lifetime of this level have been previously calculated by Hafner 
and Schwarz [31] and by Bruneau [32] . 



2 Theory 

The modified version [29] of the GRASP implementation 30J of the multiconfiguration Dirac- 
Hartree-Fock method [33] was used in the present paper. The starting point is the Dirac- 
Coulomb Hamiltonian 

H DC = «*,: • Pt + (A - i)c 2 + v t N + i/ rij , (1) 

i i>j 

where V N is the monopole part of the electron-nucleus Coulomb interaction. The wavefunction 
for a particular atomic state (<P) is obtained as the self-consistent solution of the Dirac-Fock 
equation [30] in a basis of symmetry adapted configuration state functions (#) 

NCF 

&(rPJM) = Ci${liPJM). (2) 

i 

The basis NCF was systematically enlarged [34|35] to yield increasingly accurate approxima- 
tions to the exact wavefunction. All calculations were done with the nucleus modeled as a 
variable-density sphere, where a two-parameter Fermi function [36] was employed to approxi- 
mate the charge distribution. The Breit and QED corrections were estimated with the step- wise 
procedure described in [34] . They were applied only to two transition rates, as explained in sec- 
tion [43] below. 



Will be inserted by the editor 



3 



3 Method 

The wavefunctions were obtained with the active space method in which configuration state 
functions of a particular parity and symmetry are generated by substitutions from a reference 
configuration to an active set of orbitals. The active set and the multiconfiguration expansions 
are increased systematically. The whole process is governed by convergence of the expectation 
values. The calculations were divided into two stages. Each stage was further divided into 
several consecutive steps. In the first stage, the spectroscopic and virtual spinorbitals were 
generated in relatively small multiconfiguration expansions. The spectroscopic orbitals required 
to form a reference wavefunction were obtained with a minimal configuration expansion, with 
full relaxation. Then virtual orbitals were generated in five consecutive steps. At each step 
the virtual set has been extended by one layer of virtual orbitals. A layer is defined as a 
set of virtual orbitals with different angular symmetries. In the present paper five layers of 
virtual orbitals of each of the s, p, d, f, g, h symmetries were generated. At each step the 
configuration expansions were limited to single and double substitutions from valence shells to 
all new orbitals and to all previously generated virtual layers. These were augmented by small 
subsets of dominant single and double substitutions from core and valence shells, with further 
restriction, that at most one electron may be promoted from core shells (which means, that in 
the case of a double substitution the second electron must be promoted from a valence shell). 
All configurations from earlier steps were retained, with all previously generated orbitals fixed, 
and all new orbitals made orthogonal to others of the same symmetry. The initial shapes of 
radial orbitals were obtained in Thomas-Fermi potential, and then driven to convergence with 
the self-consistency threshold set to 10 -8 . All radial orbitals were separately optimized for each 
of the nine atomic states of interest. The Optimal Level form of the variational expression |36j 
was applied in all variational calculations. 

In the second stage, the configuration-interaction calculations (i.e., with no changes to the 
radial wavefunctions) were performed, with multiconfiguration expansions tailored in such a 
way, as to capture the dominant electron correlation contributions to the expectation values. 
The valence and core- valence effects constitute the dominant electron correlation contributions 
in the oscillator strength calculations [37], therefore all single and double substitutions were 
allowed from several core shells and both valence shells (i.e., 7s 2 , 7 sip, or 7s6d, depending 
on the state) to all virtual shells, with the same restriction as above, i.e. that at most one 
electron may be promoted from core shells. The virtual set was systematically increased from 
one to five layers, until the convergence of transition rates was obtained. In a similar manner, 
several core shells were systematically opened for electron substitutions — from the outermost 
6p to 5s5p5d6s6p shells. The effects of substitutions from AsApAdAf shells were neglected. We 
estimated them separately for three states and discovered that they change the calculated values 
of transition rates by no more than a fraction of a percent. The transition rates were calculated 
with the biorthonormal technique 38 39J, which permits the application of standard Racah 
algebra, while retaining the advantage of wavefunctions separately optimized for each state. 
Experimental values of transition energies from Moore's tables [8] were used in calculations of 
transition rates. 



4 Results 

4.1 The metastable 7s7p 3 P Q state 

In principle there are three possible decay channels of the 7s7p 3 Po state. It can decay to (the 
only lower lying) ground state through (1) a blackbody radiation induced decay, (2) a 2-photon 
E1M1 transition, or through (3) a hyperfine induced transition. The first is beyond the scope of 
the present paper, since it depends on the ambient temperature. Of the other two, the former 
can be estimated through an order-of-magnitude comparison with the E1M1 ls2p 3 Po - Is 2 l So 
two-photon transition in helium-like heavy ions. We started from a recent evaluation [40] . which 
gives 3.14 x 10 9 s _1 for Ra. In order to obtain a good dependence on transition energy and 
radial matrix elements, we evaluated the E1M1 matrix elements (following Eq. (2) of Ref. [40 ) 
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for He-like Ra (using only ls2p 3 Pi as an intermediate state) and for neutral Ra (using only 
the 7s7p 3 P\). Then we integrated over photon energies with the five-point Gauss-Legendre 
formula. Finally, we evaluated the ratio of these two values to scale the He-like E1M1 rate. This 
gives an order of magnitude estimate of 9.6 x 10~ 3 s -1 . Such lifetime (« 100 s) of the 7s7p 3 Po 
state would be comparable to the lifetimes of nuclei of several radium isotopes (it would be 
significantly shorter only in comparison with the nuclei of the most stable radium isotopes 
spanning the mass range 223-229). This would be the case of spin-zero isotopes of radium. 

For the isotopes of radium with a nonzero value of nuclear spin, the hyperfine-induced 
transition must also be considered. We estimated the 7s7p 3 Pq - 7s 2 l So transition rate with a 
simple three-state model, in which the wavefunctions of the hyperfme components of the upper 
3 Pq state are described by a symmetry-adapted configuration-state-function expansion of the 
form 

\7s7p 3 P IF) HFS = c \7s7p 3 P IF) + Cl \7s7p ^IF) + c 2 \7s7p ^IF) (3) 

running over the appropriate hyperfme components IF of the l P\ and 3 P\ states. The hyperfine- 
induced transition rate (in s _1 ) may be approximately expressed as 

A(7s7p 3 P ° -> 7s 2 %) = 
2 - 0261 3 3 A 3 X 1018 \ci(7 S 7p ^WQPWr* %)+C2{7s7p *mQPP* 'Sof (4) 

where (7s7p 1 Pf|[Q^ 1) ||7s 2 %) and (7s7p 3 Pf||QP||7s 2 %) are reduced matrix elements for 
the electric dipole operator, and A is the transition wavelength (in A). The coefficients C\ and 
Ci are related to the off-diagonal mag netic dipole constants A%fi S ( 1 P 1 3 P ) and A%fi S ( 3 P 1 3 P ) 

AHFS/lp 3p \ AHFS(3p 3p N 

(see [JT] or [12] for full derivation). The sum in the expansion ([3]) should run over all excited 
states, but the sum in Eq. (Q is usually dominated by those states, for which the transition rates 
to the ground state are large and at the same time the energy denominators in Eq. ([5]) are small. 
In case of the 7s7p 3 Pq state of radium the 7s 7p 3 Pi and 7s 7p X P\ states dominate. The contribu- 
tion of the reduced matrix element (7s7p 3 P*° HQi 1 ^ ||7s 2 ^o) is 3.7 times larger than that of the 

(7s7p 1 Pf || Qx l|7s 2 ^o) matrix element. The contributions of other states are much smaller, 
due to the presence of the energy denominators in the coefficients (O and to the fact that the 
corresponding off-diagonal matrix elements are many orders of magnitude smaller (actually they 
are exactly zero, unless non-orthogonality between fully relaxed wavefunctions and correlation 
effects arc not neglected), and they can be safely ignored at present level of the overall accuracy. 
The rate of the hyperfine-induced transition is isotope-dependent, i.e. it depends on the nuclear 
spin and on the nuclear magnetic moment. As in our previous paper [27], the 2 ||Ra isotope was 
chosen to set the nuclear parameters (the transition rate A( 3 Po — So) may be readily recalcu- 
lated for other radium isotopes, for which nuclear spins and magnetic moments are known). The 
nuclear spin of 2 ||Ra is I — 3/2, and the nuclear magnetic dipole moment fj, = 0.2705(19)/iat 
was taken from the paper of Arnold et al [13 . The electric dipole transition rates AQ-P\ — So) 
and A( 3 Pi — ^o) from tableQ](in Babushkin gauge) were used. The off-diagonal magnetic dipole 
constants A^^P^Pq) = 540 MHz, and A^ s ( 3 Pi, 3 P ) = 1172 MHz, were evaluated with 
the use of the same wavefunctions, and in the same approximation, as described in section ([3|) 
above. Together they yield hyperfine-induced transition rate ^4( 3 Po — £>o) = 0.0210 s _1 in case 
of constructive interference, and A( 3 Po — &o) = 0.0070 s _1 in case of destructive interference. 

The above evaluations of the E1M1 and hyperfine-induced rates were performed indepen- 
dently. As it turned out, the two contributions are of the same size, therefore they have to 
be treated simultaneously. To this end we employed the effective Hamiltonian method [43]. It 
requires also an evaluation of the A^ s ( 3 F'i, 1 Pi) matrix element (851 MHz). The method is 
valid beyond the limits of the perturbation method exposed above, in particular it does not 
require that the energy separation between levels is large compared to level widths. It yields 
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A ( 3 Po — iSo) = 0.02935 s 1 (if the two-photon transition is switched off the effective Hamilto- 
man method yields A ( 3 P - &o) = 0.0197 s" 1 ). In the process of this calculation, we found out 
by comparison with the results of mdfgme code [44j that non-orthogonality between spinorbitals 
plays very important role in the evaluation of non-diagonal hyperfine matrix elements (neglect- 
ing it may even lead to a sign change). We employed a new code, developed by one of us (P. J.) to 
evaluate correlated off-diagonal hyperfine matrix elements from GRASP wavefunctions, taking 
into account the effect of non-orthogonality between spinorbitals. 

The accuracy is limited by the electric dipole matrix elements (7 's7p 3 P 1 °||Q^ 1 ' ) ||7s 2 ^So) and 

(7s7p ipfHQ^ ||7s 2 ^o) as well as by the off-diagonal magnetic dipole constants Af£ s (^Px^Pi), 

Af f f s ( 3 Pi, 3 P ) and A^^P^Po). In case of the matrix element (7s7p 3 P J 1 °||Q^ 1) ||7s 2 %} we 
may assume the 5% relative accuracy for the 3 Pi — transition rate from experiment [5J- The 
accuracy of the A^Pi — 1Sq) rate is more difficult to estimate (see ref. [IS]). This is the strongest, 
'allowed' transition in the radium spectrum, but (as mentioned above) the contribution of its 

matrix element (7s7p l P°\\Q,i \\7s 2 ^So) to the total value of the calculated rate of the hyperfine- 

induced transition 3 Po— ^Sq is 3.7 times smaller than that of the (7s7p ^"HQ^I^s 2 ^So) matrix 
element, so even a relatively large error bar would be quenched. We may very conservatively 
take the entire difference between the Babushkin and Coulomb gauge final values from table 1 
in reference [5SJ as the error limit, obtaining a relative accuracy 25% for the A^Pi — x Sq) 
transition probability. The accuracy of the diagonal magnetic dipole constants was estimated 
to be 6% (see [26 ). Since an off-diagonal constant depends on both states rather than on one, 
we, again quite conservatively, doubled the 'diagonal' error limit estimate and assumed 12% 
as the contributions of the off-diagonal hyperfine constants to the error bar. Eventually, the 
above procedure yields the lifetime of the metastable 3 Po state t = 34(15)s of 2 gf Ra isotope. 
The lifetime is based on both the 2-photon and hyperfine-induced channels. The error bar does 
not include the 2-photon contribution. 



4.2 The 7s7p 3 Pi state 



Figure Q] presents the transition probability A( 3 Pi — 1 Sq) as a function of the size of the mul- 
ticonfiguration expansion. The transition rates were calculated in Babushkin and Coulomb 
gauges, which in the non-relativistic formulation correspond to length and velocity form of the 
transition integral, respectively. Both curves were obtained in the 'core-valence' approximation 
described above. The resulting Babushkin and Coulomb values arc compared with the experi- 
mental result obtained by the Argonne group [5] and with three available theoretical values (the 
value A=A ■ 10 6 s _1 obtained by Bruneau [32] did not fit within the vertical scale of figured]). 
It is clearly seen in the figure Q] that the core- valence correlation effects are saturated in the 
framework of the single and restricted double substitutions and five layers of virtual orbitals, 
as described in section ([3]) above. The remaining difference between the final Babushkin and 
Coulomb gauge values may be attributed to the omitted core-core effects. We have made an 
attempt to estimate the contribution of the core-core correlation to 3 P/2 — Sq and x Pi — ^So 
transition rates [25j , but for other transitions the gauge differences and comparisons with data 
obtained by other authors, where available, are the only indications of the accuracy of our 
calculated rates. Although the gauge difference must not be treated as the error bar per se, it is 
a useful indicator in partially saturated multiconfiguration calculations of transition rates. The 
values obtained in the Babushkin gauge are weighted toward the outer parts of the electronic 
wavefunctions, while the Coulomb gauge values weight more inner parts, where the core-core 
effects arise. Therefore partially saturated expansions often produce Babushkin and Coulomb 
gauge values converging toward different limits, as in figure [T] with the difference arising from 
the omitted core-core effects. This is also the reason why Babushkin gauge results are usu- 
ally treated as more reliable, which seems to be confirmed by the agreement of the Babushkin 
gauge transition rate A( 3 Pi — So) with experiment (although agreement this close is very likely 
accidental; and good agreement for one level is not enough to justify a more general rule). 
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Fig. 1. Transition probability A( 3 Pi — 1 So) in Babushkin (upper curve with circles) and Coulomb 
(lower curve with squares) gauges, as a function of the multiconfiguration expansion (in core-valence 
approximation — see text for details). The lone dot with error bars (to the right from the end of 
the Babushkin curve) represents the experimental result from reference [5]- The three stars at the far 
right represent the theoretical data obtained by other authors and are denoted by publication year in 
brackets: [1978] — reference [3T]; [2000] — reference g5]; [2006] — reference [25]. 



4.3 The 7s6d 3 D 2 and 7s7p x Pi states 

These two levels are distinguished because they were the subject of a separate paper [25]. In the 
present paper we duly quote the data obtained in core-valence approximation, and the reader 
is referred to the above mentioned paper for further discussion. 



4.4 The remaining 7s7p 3 P 2 , 7s6d 3 D lt 7s6d 3 D 3 , and 7s6d 1 D 2 states 

The remaining 7s 7p 3 P 2 , 7s6d 3 Di, 7s6d 3 D%, and 7s6d 1 D 2 states, together with the states 
discussed in sections |4.1[ 14.21 and 14.31 above, constitute full set of states arising from the three 
lowest electronic configurations of radium. 

Some of these levels should also be considered metastable, not only 7s7p 3 Pq- The 3 D^ — 3 D 2 
transition is the strongest 'direct' decay channel for the 3 D% state, but it is in fact very weak 
(comparable to that of the 3 Pq — transition). The rates of other possible 'direct' decay 
channels would be still smaller (e.g. M2 transition 3 Ds — 3 P\ rate is comparable to that of the 
3 D 2 — 3 Pq transition) , therefore the multiphoton or hyperfine-induced transitions may also play 
a role in the 3 D^ lifetime. Similar considerations may in principle apply to the 3 P 2 , 1 D 2 , and 
3 D 2 states (we did not pursue this issue, though). 



5 Summary and conclusions 

Table Q] presents calculated transition rates, and the lifetimes are presented in tabled] Coulomb 
gauge values of electric multipole transitions from the table Q] were not used to obtain the 
lifetimes. With the exception of the 3 Po level, the transition rates are given with 4 significant 
digits and the lifetimes with 3 significant digits, but that not necessarily reflects their accuracy. 
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Table 1. Calculated transition rates between nine lowest levels of radium [s -1 ]. Transition multipo- 
larities are denoted by El, E2, Ml, M2; the Babushkin and Coulomb gauge values by (B) and (C), 
respectively. HFS means hyperfine-induced. Numbers in brackets represent powers of 10. 





transition 


This work 


Ref. [28] 


Expt. [5] 


3 Po- 


So 


HF^+P1 A/T1 


2.935[-2] 






3 P1- 




E1(B) 


2. 374 [6] 
1.873[6] 


2. 760 [6] 


2.37(12) [6] 


3 P1- 




Ml 


1.334[-2] 






3 P1- 


3 Dl 


v J 

(C) 


8.794[1] 
4.025[3] 


9.850[1] 




3 P1- 




TjI IDI 

(C) 


1.775[-3] 
1.607[2] 


1.572[-3] 




3 P 2 - 


3 Pn 
11) 


F9fR^ 

JZjZ l J3 1 

(C) 


1.185[-2] 
1.243[-2] 






3 P 2 - 




(C) 


4.310[3] 
9.722[3] 


4.897[3] 




3 P 2 - 


3n, 
U 2 


Hi ± I _D 1 


4.602[4] 


5. 204 [4] 






3 n, 


(C) 


1.109[5] 






3 P 2 - 


\ ) 


1.044[5] 
4.201 [5] 


1.234[5] 




1 Fl- 


1 Q 




1.793 [8] 


1.805 [8] 






3 Dl 


(M 


1.795 [8] 






X P1- 


E1(B) 


3.282[4] 
5.222[4] 


4.195[4] 




L P1- 


3 D 2 


E1(B) 


9. 793 [4] 
1.441 [5] 


2. 646 [4] 




L P1- 


X D 2 


E1(B) 
(C) 


3.241 [5] 
5.875[5] 


3.194[5] 




3 Di- 


3 Po 


E1(B) 
(C) 


1.390 [3] 
7.940[3] 


1.529 [3] 




3 D 2 - 


% 


E2(B) 
(C) 


2.524[-l] 
1.630[-1] 


3.032[-l] 




3 D 2 ~ 


3 Po 


M2 


3.021 [-13] 






3 D 2 ~ 


3 d 


Ml 


5.082[-4] 






3 D S - 


3 D 2 


Ml 


6.352[-3] 






l D 2 - 


'So 


E2(B) 
(C) 


2.710[1] 
2.271[1] 






1 D 2 - 


3 Pi 


E1(B) 
(C) 


6. 960 [2] 
1.224[3] 


7.722[3] 




1 D 2 - 


3 P 2 


E1(B) 
(C) 


5.930[0] 
7.535[0] 


7.973[0] 





As discussed in section (|4.2[) above, the accuracy of our results is difficult to estimate, but it is 
probably much worse than 3 or 4 digits, particularly for very weak transitions. 

The results of the present calculations may be considered as fully converged in core- valence 
approximation (convergence has indeed been observed for all transitions, similarly to that shown 
in figure Q}, with core-core effects omitted. 

With the exception of the l P\ — 3 D2 and X D% — 3 Pi transitions (we cannot offer any plau- 
sible explanation for these discrepancies), our Babushkin gauge values are in reasonably good 
agreement with the results of Dzuba and Ginges |28j . An interesting feature is the large discrep- 
ancy between the results obtained in Babushkin and Coulomb gauges for transitions connecting 
closely- lying triplet P and triplet D states. The B/C ratios (i.e. the ratios of Babushkin versus 
Coulomb gauge results) turned out to be much closer to unity for those transitions when the 
experimental energies in the transition operator were replaced by theoretical ones. This obser- 
vation, together with inconsistent results of the calculations of energies of excited states, had led 
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Table 2. Calculated lifetimes of eight lowest excited states of radium, compared with data from other 
authors. 



state This work Ref. 



Ref. \M Ref. ED 



Ref. 



Expt. [5] 



3 Po 


34(15) s a 








421 ns 


362 ns 


505 ns 




6.46 us 


5.55 us 


5.2 us 




5.56 ns 


5.53 ns 


5.5 ns 


3 Dx 


719 ns 


654 us 


617 us 


3 D 2 


3.95 s 


3.3 s 


15 s 


3 D 3 


157 s b 






'D 2 


1.37 ms 


0.129 ms c 


38 ms 



420 ns 250 ns 422(20) ns 



(a) based on E1M1 and hfs-induced decay channels in the ^ggRa isotope 

(b) based on 3 L>3 — 3 Z>2 magnetic dipole transition only; other possible decay channels neglected 

(c) calculated from transition probabilities quoted after Table VI in Ref. [28] 



us [35] to suggest an experimental verification of radium data in Moore's tables [Sj. At that time 
we were not aware of the papers by Eliav et al [53] and Landau et al [21] , where the excitation 
energies of radium were calculated in the coupled-cluster approximation. More recently these 
calculations were independently confirmed within the framework of the CI+MBPT theory 28] . 
In both cases good agreement with experiment had been achieved. There are two principal 
differences between the methods and approximations used in the abovementioned three papers 
with respect to the methods and approximations used in the present paper. The transition rates 
were the primary targets of the present calculations, not the transition energies. Therefore we 
optimised the electronic wavefunctions separately for each of the nine states of interest. The 
core- valence correlation effects were fully accounted for. The core-core correlation effects, which 
are less important in calculations of transition rates and hyperfine structures, were treated in a 
very crude approximation in the cases of So, 3 I?2 and 1 Pi states, and were neglected for other 
states. The calculations of energy level differences require well balanced orbital sets and exten- 
sive multiconfiguration expansions. The results are usually in better agreement with experiment 
if a common set of orbitals is used for both states. If the orbital sets are separately optimised, 
the transition energy is obtained as a pure difference of the total energies of the two states 
of interest. They are both several orders of magnitude (five orders in case of radium) larger 
than the transition energy itself, therefore our calculated transition energy values may be less 
accurate than the calculated transition rates and hyperfine structures. Any further refinement 
of the present calculations would require computer resources, which are currently unavailable. 
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